Vibrating Superconducting Island in a Josephson Junction 
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We consider a combined nanomechanical-supercondcuting device that allows the Cooper pair 
tunneling to interfere with the mechanical motion of the middle superconducting island. Coupling 
of mechanical oscillations of a superconducting island between two superconducting leads to the 
electronic tunneling generates a supercurrent that is modulated by the oscillatory motion of the 
island. This coupling produces alternating finite and vanishing supercurrent as function of the 
superconducting phases. Current peaks are sensitive to the superconducting phase shifts relative to 
each other. The proposed device may be used to study the nanoelectromechanical coupling in case 
of superconducting electronics. 
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Vibrational modes (vibrons) and spins possess dy- 
namical degrees of freedom and they have a large im- 
pact on electron dynamics. Peaks and dips in the dif- 
ferential conductance of molecular electronics devices 
[DI3iaffl0ianEiai]lini]2II3mmay indi- 
cate strong effects from electron-vibron coupling. Steps 
in the differential conductance have been observed, both 
experimentally and theoretically, in STM based inelastic 
tunneling spectroscopy (IETS) around local vibrational 
mode on surfaces [TSl US] • Effects from local vibrational 
modes on the conductance in molecular quantum dots 
and single electron transistor have also been investigated 

[m mi na uni nn na na mi nu- These studies 

may have implications on charge-based quantum infor- 
mation technology. 

More recently the field has been moving towards in- 
corporating superconducting electronics into nanoelec- 
tromechanical devices, e.g. see recent work [28l|29]. One 
example of superconducting electronics in combination 
with nanomechanical setup is a Cooper pair shuttle. The 
setup proposed originally [SU] was aimed at a steady state 
description of the mechanical assisted Cooper pair tun- 
neling. We extend the analysis of Cooper pair shuttle to 
consider possible resonances between mechanical motion 
and ac Josephson effect. We here address dynamical as- 
pects of this device where we focus on time domain in 
the presence of ac signatures in the electric current. 

In this Letter we study Josephson tunneling in a sys- 
tem with a mechanically moving superconducting island 
between two superconducting leads. Mechanical oscilla- 
tions of the island couples to the electronic tunneling. 
This coupling gives rise to an oscillator modulation of 
the Josephson current, such that the Fourier spectrum 
in general exhibits current peaks at loj ± loq, loj ± 2ujq, 
and 2luj, where loj =eV and ujq are the Josephson and 
renormalized oscillator frequencies, respectively. In addi- 
tion, the zero bias supercurrent peaks at the frequencies 



0, ±ujq, and ±2£o, however, the peaks at and ±2iI>o, 
or ±iI>o vanish for certain combinations of the supercon- 
ducting phases in the three components of the system. 

The physical system we are consider consists of a su- 
perconducting island between two superconducting leads, 
where the island is exposed to mechanical vibrations 
modeled with Hook's law force constant k c , see Fig. [I] A 
bias voltage is applied across the junction. The Hamil- 
tonian for the system is given by 

H = Hl+Hi + H r + H t . (1) 

Here the leads (L, R) and the island (I) 
are described by the BCS Hamiltonians as 

^■L(I,R) = Y J p(k,q)<j e p{k,q) C p (k,q)a C P{k,q)<? + 

E p ( k ,g)[ A HI,R) C l{k, q n C -p{-k,- q )l + H - C -l- We den0tc thc 

creation (annihilation) operators by c , k ^ a (cp(k.q)a), 
where the subscript p (fc, q) is momenta in the leads 
and island, whereas a =tjj, is a spin index. Finally, 
£ p(k,q) an d ^l(i.r) are the single-electron energies and 




FIG. 1: Schematic view of the mechanically and electronically 
coupled superconducting island (SC I) to the superconduct- 
ing leads (SC L and SC R). The cantilever superconducting 
island is modeled as a harmonic oscillator with spring con- 
stant k c and mass m c , placed between the infinitely massive 
superconducting leads. The device is biased with voltage V . 
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pair potential (gap function), respectively. Without 
loss of generality, we assume that the superconduc- 
tors are of a conventional spin-singlet s-wave pairing 
symmetry, and consider the Josephson tunneling at 
zero temperature. The last assumption implies that 
the temperature is low compared to relevant supercon- 
ducting gaps, e.g. T <C A L ( IR y To avoid thermal 
damping of an oscillator we also require low enough 
T ^ 10 — 100 mK. The last term in Eq. (JIJ describes 
the tunneling between the leads and the island, i.e. 
n T = I2pkA T pkcU C k<? + H - C A +J2 q k*\ T ikc\ a c krT + H.c], 
where the tunneling matrix elements T p ^ k transfers 
electrons through insulating barriers between the leads 
and the island. The local vibrational mode of the island 
is in the linear coupling regime given by 



T pk = T$(l + a L u), T qk = T$(l + a R u), (2) 



frozen (u = 0), given by 



E n E k \uj x + E n + E k 



— E n — E k J 



(4) 



where A x and Aj is the superconducting gap in the 
lead x an d island, respectively, whereas the quasi- 
particle energies E n = y/ (e n — ^ x ) 2 + |A X | 2 and E k — 
\J (e k — /j,i) 2 + |A/| 2 . The second contribution to the 
Josephson current in Eq. ^ has the amplitude 



E n E k \(uj x + E n + E k ) 2 



{u x -E n - E k y 



(5) 



where ocuja describes the coupling between the tunnel- 
ing electrons and the vibrational mode corresponding to 
the left (right) tunnel junction. The quantity u is the 
displacement operator for the oscillator. The tunnel- 
ing matrix clement T pk is exponential in displacement u, 
thus, the assumed linear coupling is a good approxima- 
tion for small u. This allows evaluation of aL.R in terms 
of the tunneling matrix elements and their distance de- 
pendence. We assume here a very general equilibrium 
geometry with no particular symmetry being required. 
The equilibrium point (u = 0) of the mechanical oscilla- 
tor placed within the junction could be placed anywhere 
in between the leads. The energy associated with the 
vibrational mode, uiq = ^k c /m c ~ 10 _1 — 10~ 6 eV, is 
much smaller than the typical electronic energy on the 
order of 1 eV, the mechanical oscillations are very slow 
compared to the time scale of the electronic processes. 
This allows us to apply the Born-Oppcnhicmcr approxi- 
mation to treat the electronic degrees of freedom as if the 
local oscillator is static at every instantaneous location. 

The current is derived using standard methods. For a 
given bias voltage eV = Hl — Mfl; where /i x , \ = L, R is 
the chemical potential of the left (L) and right (R) lead, 
the Josephson current between the lead \ an d the island 
is given by [S7J (setting H = 1) 



J|(t) = J*(u; x )[l + a x u} 2 sm(u J }t + <f> x ) 



-rg(« x )[i + . 



Y x ii cos (ujjt + cf> x ) , (3) 



in applying the local approximation u(t') ~ u(t) + (t' — 
t)u(t). Here we have introduced the phase difference 4> x 
between the lead x an d the island, and the Josephson 
frequency = 2(fi x — fij). In Eq. Jg is the am- 
plitude of the Josephson current when the cantilever is 



Using that Ig = —Ig for stationary bias voltages we 
write the total Josephson current Ig = Ig = (Ig — Ig )/2. 

The Hamiltonian Hj for the Josephson energy [27J is 
derived by requiring that the derivative of TLj with re- 
spect to 4> x yields the supercurrent given in Eq. (j3). We 
find that 



\ ( y E}(l + a x u) 2 [l-cos(u J x J t + 



X=L,R 
1 

2~e 



Ig(l + a x u)a x ii sin(ujjt + <p x ) 



(6) 



where E* = Jg/(2e). This effective Hamiltonian cap- 
tures the back action effect and reflects how the electronic 
degrees of freedom influence the mechanical oscillator. 
The classical motion of the shuttling island is described 
by the Hamiltonian 7i osc = p 2 / (2m)+k c u 2 /2+Hj , which 
results in the classical equation of motion 

m c u+ [75 (t) + 7jv]u + k c u = F(t). (7) 

Here, the driving force 



F(t) 



£ E*a x (l + a x u) l-(l + ^ 



X=L,R 

x cos(u)jt + <f> x ) 



, .x-px 
AeE X j 



(8) 



whereas the time-dependent damping factor 75 (t) + 77V 
contains the superconducting contribution 

7s(*) = -^ E rga x sin(^ + ^ x ). (9) 

X=L,R 

due to damping energy in and out of the superconducting 
carriers, and the external damping 7jv- We assume that 
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the damping into non-superconducting degrees of free- 
dom is suppressed by the factor exp (— Al^/T), as one 
would need to damp energy into normal excitations. 

In order to make next steps analytically, we assume the 
system to be in the weak coupling limit. Moreover, since 
<C Jg the main physics is captured by neglecting the 
terms in the damping and driving force proportional to 
a x T^ and a 2 . The motion of the island is then given by 



u(t) = u (i sm(u t + 8 a )e" lNt/2m " 



x{cos (u)jt + fa) + B x sin (upjt + fa)}], 



where 



1 - (u^/wo) 5 



[1 - (u^M,) 2 ] 2 

jNOJj 1 



(-y N uyk c y 



i - K/w ) s 



(10) 

(11) 

(12) 



The first term describes the unperturbed motion of 
the island around its equilibrium position, where Uq 
and are to be determined by the initial conditions, 
whereas e± = — 7at/2to ± iuj$ with cigenfrequcncies 
ujq = \/ lUq — [7„/2m] 2 of the mechanical oscillations. 



The other terms arise due to the coupling between the 
mechanical and electronic degrees of freedom. Using this 
expression for the mechanical motion of the island, we 
obtain the Josephson current approximately as 

i|(t) = J|[l + 2a x sin(£ t + S )e'^ Nt / 2m - 



+a 2 x sin 2 {Q t + 5 )e-^ t l m ^2a x £ 

X'=L,R x ' 

- A x ,{cos(u;X't + <t> x ,) + B x , sm^'t + fa)})] 



x sin(w^f + fa), 



(13) 



where a x — uqol x and K x — k c UQ/E x . 

We now consider a few limiting cases in the symmetri- 
cally biased system, such that ojl = — = lo,j (= eV), 
and we assume that Tpn = Tq„ and |Ax,| = \Ar\. Then, 

jf 



jR 



Js, hence, we set Ej 



K. 



Ej and K x 

Moreover, Al = Ar = A, and Bl = —Br = B, and we 
also assume that ckr, = —a^ = —a. Using Eq. (13 1 we 



find that the total current under those circumstances can 
be written as 



2ig(*) 
Js 



[1 + a 2 sm 2 (£j t + S )e lNt/m "][sin(ujjt + fa) + sm(uj, ; t - fa)] + 2Ssin(u> * + 5 )[sm(ujjt + fa) 

- sm(ojjt - (t} R )]e- lNt/2m - + 5 2 yl[sin 2{ujt + fa) + sin 2(ujjt - fa) - 2 sm(2ujjt + fa - fa) 

+ B{2 ~ cos 2(uj. ; t + fa)- cos 2(ujt -fa>) + 2 cos(2uijt + fa - fa) - 2 cos(<j) L + fa)}]/K (14) 



In the case where the phases 
rent thus reduces to 

2I S (t) 
Js 



2 l+a 2 sin 2 (u} t+5 )e 



the Josephson cur- 



^ f/ "'j smujt. (15) 



In the undamped, or very weakly damped case, this ex- 
pression predicts that there will be a dc component to the 
Josephson current whenever the bias voltage matches the 
Shapiro step value 2w$. On the other hand, there is no 
dc component to the current at the value (Do as is pre- 
dicted for a single tunnel junction in Ref. [27 . The 
motion of the island between the two leads, clearly, has a 
compensating effect in the sense that the coupling of the 
mechanical oscillator and the tunneling electrons cancel 
their respective side peaks at loj±loq- These observations 
open up possibilities to excite cantilever motion by gener- 
ating the Josephson current with appropriate frequencies. 



the absence of Josephson current at ujj ± lvq and 2u> j, 
while there is a finite Josephson current at lo,j ± 2<I>o , are 
clearly illustrated in Fig. 2 These plots show the Fourier 
transform of the Josephson current log 10 \2Is(u>) / Js\ 2 , in 
the undamped case (7^ = 0), as function of the phase 
fa and frequency u, for the fixed phase 4>r — and two 
different frequencies uiq. 

More interesting is the zero bias Josephson current 
found from Eq. ( 14 1 . Assuming that the phases <f>L — 
(bp = 6 we obtain 



21 s {t) 
Js 



4Ssin(w o t + (So)sin0e~ 7JV * /2mc (16) 



The features predicted from Eq. (14 1 at fa = 0> l na ^ ^ s 



Hence, there is an ac component which is modulated by 
the eigenfrequency ujq of the moving island. The dashed 
lines in Fig. |3]signify equal phases fa- = (f>, at which only 
the ac component which is modulated by luq is finite. 
In the case of opposite phases, 4>l = —<j)R — 4>, we find 



4 



9- 




FIG. 2: (Color online) Fourier transform of the Josephson 
current log 10 |27s(a;)/Js| 2 in the symmetrically biased and 
undamped system as a function of the frequency uj and the 
phase 4>l, for a fixed phase <f>n = 0. Here we have taken a = 
0.1, K = 0.6, and the oscillating frequencies a) luq — OAloj 
and b) ujq — lAuo.j. 



the Josephson current 



Js 



= 2 



(\ + a 2 sm 2 (u t + 5o)e-" tNt/m ^j sin cj> (17) 



In contrast to the previous case, in this case there is a dc 
component depending on the phase (f>. Moreover, oppo- 
site phases give rise to an ac component which is mod- 
ulated by twice the eigenfrequency (2loq) of the moving 
island. Opposite phases are signified by solid lines in 
Fig. [3j showing a finite dc component to the Josephson 
current, along with finite ac components at ±2<I>o. 




FIG. 3: (Color online) Zero bias Josephson current 
log 10 \ 2Is{co) I Js\ 2 in the symmetrically biased and undamped 
system as a function of the frequency ui and the phase </>£, for 
a fixed phase 4>R = tt/2. Here we have taken ujq = 1, while 
other parameters are as Fig. [2] 



For an experiment of this type, taking |A| ~ 10 meV 
as relevant to MgB2 [31] provides a sufficiently small me- 
chanical damping. By acquiring a vibrational frequency 
of wo/27r ~ 1 GHz [32], it should be possible to tune the 
Josephson frequency such that ujq/luj > 0.1. A coupling 
strength a/u} ~ 10 _1 — 10~ 3 [33J would be sufficient to 
enable read-out of our predictions. 
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